Verifying multi-partite mode entanglement of W states 
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We construct a method for verifying mode entanglement of TV-mode W states. The ideal W state contains 
exactly one excitation symmetrically shared between TV modes, but our method takes the existence of higher 
numbers of excitations into account, as well as the vacuum state and other deviations from the ideal state. 
Moreover, our method distinguishes between full TV-party entanglement and states with A/-party entanglement 
with M < TV, including mixtures of the latter. We specialize to the case TV = 4 for illustrative purposes. In the 
optical case, where excitations are photons, our method can be implemented using linear optics. 
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I. INTRODUCTION 

Detecting and classifying entanglement is an important 
challenge in the field of quantum information science. One 
problem is of a theoretical nature, to decide whether a given 
density matrix p of multiple quantum systems is entangled or 
separable. Even for bipartite systems, this is a hard problem 
for which no efficient general solution is known for higher- 
dimensional Hilbert spaces, although a simple test based on 
the negativity of the partial transpose of the density matrix 
leads to a sufficient criterion for entanglement |l j. If p is en- 
tangled, the next issue is how to classify the type of entangle- 
ment. For more than two subsystems, the full classification of 
all entanglement classes is as yet an unsolved problem (see, 
e.g., Refs. 

In an experiment the practical task of detecting entangle- 
ment is even harder. If one would perform a full tomographic 
measurement, then in the limit of infinitely many data one 
would end up with an arbitrarily accurate estimate of a den- 
sity matrix p, and thus reduce the experimental problem to 
the theoretical problem mentioned above. In all other cases, 
one needs different tests that make use of less than full knowl- 
edge of the density matrix. The main practical disadvantage 
of full tomography is the rapidly growing number (with the 
number of quantum systems and with the dimension of the 
Hilbert spaces involved) of measurements needed to find all 
elements of a density matrix. The other challenge is obtain- 
ing a physical density matrix from a finite set of measured 
data El E| . 

Thus there is an evergrowing demand for simpler experi- 
mental tests revealing entanglement. Fortunately, for bipartite 
systems there exists a handful of different experimental tech- 
niques for entanglement detection [6J O . Here we will focus 
on a particular type of multi-partite entangled states (namely, 
W states [2|) that can be produced in systems with variable 
numbers of excitations. We think in particular of experiments 
with atomic ensembles (see, e.g. HI) based on the DLCZ pro- 
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tocol (5] in which information can be stored in the number 
of atomic excitations of each ensemble, as well as of experi- 
ments on photonic systems (see, e.g. ifTOl ), where the number 
of photons in a given mode can be used as a quantum vari- 
able. In the following we will use the words "excitation" and 
"photon" interchangeably. 

We define the state \W) as a mode-entangled analogue of 
standard TV-partite W states of qubits. It is a pure state where 
a single excitation is shared symmetrically between TV modes 
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where |0) denotes a state of a mode with no excitations and 1 1) 
is a state with a single excitation. The subscripts i = 1 ... TV 
refer to modes that are in spatially distinct locations, so that 
the concepts of "local operations" and hence entanglement are 
unambiguously defined ifTTI . 

We solve the problem of detecting the entanglement of a W 
state (and it noisy cousins) in two steps. In Section II we will 
focus on detecting and classifying entanglement within the 
subspace of a fixed total number of excitations (in all modes 
together), namely one. In Section III we complete the analy- 
sis by including the remaining parts of Hilbert space, the sub- 
space with no excitations and the subspace with more than a 
single excitation in total. Including both subspaces is crucial 
in the analysis. Earlier detection schemes [12| for W states in 
the context of photons were incomplete due to the neglect of 
states with multiple photons. Moreover, we will discuss how 
to include imperfections such as losses, most releveant for the 
actual implementation of our method lfl"3l . 

II. GENUINE TV -MODE ONE-PHOTON ENTANGLEMENT 

TV parties can be entangled in many different ways. In some 
papers "genuine" TV-party entangled states include states that 
are mixtures of M -party entangled states with M < TV, as 
long as such mixtures are not biseparable along any particular 
splitting of the TV parties into two groups (for instance, O). 
Here, however, we will classify such mixtures as A/ -party en- 
tangled states, and the name "genuine TV-party entanglement" 
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in our case is reserved for states that can only be written as a 
mixture of pure states that all possess TV-party entanglement. 
Thus our criterion for genuine TV -party entanglement is more 
severe. 

Recently it has been suggested that uncertainty rela- 
tions can be used as an entanglement criterion for finite- 
dimensional systems ( ifTSl : see also 1 16 1 for an experimental 
implementation using local observables on two qubits). The 
uncertainty principle sets up a fundamental limit on how accu- 
rately observables of a quantum system can be simultaneously 
determined. For instance, if {Mi}, i = 1 ... K is some set of 
observables, then the measurement uncertainty in a given state 
p is given by the sum of variances of all observables Mj, i.e., 
J2i 5Mi{p) 2 . This sum equals zero if and only if the state for 
which measurements of all Mj are performed is a simultane- 
ous eigenstate of all M,. If there is no such state (when the 
observables are not all mutually commuting) then there exists 
a positive number C such that 
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(2) 
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Hofmann and Takeuchi pointed out that the existence of the 
lower uncertainty bound C can be employed as a separability 
criterion 1131 . Indeed, if for some fixed set of observables an 
inequality of the form (|2j holds for all separable states then its 
violation is a signature of entanglement. 

The uncertainty bound has another obvious but important 
property. Namely, one can never decrease the average uncer- 
tainty by mixing different states. In other words, for any state 
P = J2mPmPm, an d any observable A the following inequal- 
ity holds 



SA( P ) 2 > Y,PmSA( Pm ) 2 



(3) 



The proof is rather straightforward and can be found in JT3J. 

With the uncertainty criterion at hand we still have some 
flexibility over the type of observables to choose. In princi- 
ple, all observables can be divided into two groups - local and 
nonlocal. Whereas local observables can be measured sepa- 
rately for each and every party and therefore tend to be easier 
to access in an experiment, they often cannot reliably detect 
genuine multipartite entanglement. Nonlocal observables, on 
the other hand, require a simultaneous nonlocal measurement 
of several parties at a time, which often is experimentally chal- 
lenging. Here we show how experimentally accessible nonlo- 
cal observables can be constructed to unambiguously detect 
genuine multipartite entanglement of the W-type. 

The basic idea behind the construction of nonlocal observ- 
ables is to choose them as projectors onto a basis of TV-partite 
entangled states. Simultaneous eigenstates of these projec- 
tors are necessarily entangled states, and the variance in the 
projectors is minimized for TV-party entangled states. A suffi- 
ciently small variance is then a sufficient criterion for genuine 
TV-party entanglement. In order to illustrate this idea we will 
consider a system of four modes sharing a single photon. The 
problem at hand is then to find a set of nonlocal observables 
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FIG. 1: Beam-splitter setup to project onto the four W states 
4 input modes are converted into 4 output modes by 4 50/50 lossless 
beamsplitters, numbered 1-4. From the count statistics of (ideal) 
detectors placed at the 4 output modes one obtains the quantity A(p) 
defined in Eq. ([8j. Losses, asymmetries in the beamsplitters, and 
non-ideal detectors are discussed in SectionlTv] 



which allows to separate all four-mode separable and bisep- 
arable states from the genuinely four-mode entangled states 
such as the W state of Eq. ([TJ. We note that the general con- 
struction for an arbitrary TV can be done in a similar fashion. 
Moreover, the nonlocal observables for single photons we use 
can be measured just using linear optics (beam splitters) and 
non-number resolving photodetectors. 

The Hilbert space of a system of four modes sharing ex- 
actly one photon is spanned by four basis product vectors 
{ 1 1000), |0100), |0010), |0001)}. This basis can always be ro- 
tated to a basis constituted by four W-like states, 
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The mode transformation from the four product states to these 
four W-like states can be easily decomposed in terms of uni- 
tary operations that can be implemented with beamsplitters 
and phaseshifters (see Figure [TJ. 

The next step is to choose four projectors onto the basis 
Eqs. QTJ as nonlocal observables Mi = \Wi){Wi\,i = 1,4. 
Clearly, the only simultaneous eigenstates of all four operators 
Mi are the four states | Wi). The total variance of all Mi's van- 
ishes for any one of the states | Wi). In contrast, for product 
states the total variance is bounded from below, since there ex- 
ists no simultaneous product eigenstate of all the Mj. There- 
fore, we can write down an uncertainty-based entanglement 
criterion using nonlocal observables for any state pi within 
the subspace of a single excitation, in terms of the sum of 
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variances of Mi, 

4 

A( P1 ) = $3li(pi[|Wi)<Wi|] a )-rrr(pi|Wi)(Wi)] 2 

i=l 
4 

= [Wi/'iiwi) - W*|pi|Wi> a ] 
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4 

= l-^^lpilWi) 2 , (8) 

i=l 

where the subscript 1 is there to remind us the state contains 
exactly 1 excitation. 

To find the lower bound on A for unentangled states it is 
sufficient to consider pure states thanks to Eq. d3j. For a pure 
state pi — | a) (a | we have 
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The next step is to find the minimum of A(pi) by maximizing 
J2i I (Wt|a) | 4 over all separable states |a) containing a single 
excitation. There are three types of pure four-mode states that 
are not four-mode entangled GTl 

1. The fully separable pure states, which are products of 
four single-mode states. There are only four such states 
within the subspace of interest, namely 1 1000), |0100), 
... |0001). 

2. Biseparable states with at most two-mode entangle- 
ment. Here two modes must be in the vacuum state, 
and the most general pure state in this class is of the 
form 1 00) ® (a|01) + 6| 10) ), or similar states resulting 
from permuting the different modes. 

3. Biseparable states with at most three-mode entangle- 
ment. Here at least one mode is in the vacuum state, 
and the most general pure state, up to permutations of 
modes, is of the form |0) ® (a|001) + 6|010) + c|100)). 

Given the most general pure state within each class it is 
straightforward to calculate the three corresponding minimum 
values of A(pi), and the results are depicted in Fig. [2] For 
example, for any pure fully separable state \a) the overlap 
KWila)! 2 = 1/4 for any i, and so A(pi) =3/4. For general 
mixtures of fully separable states, this number gives the best 
possible lower bound on the variance. We note that the nu- 
merical results from the next Section will confirm the results 
of Fig. [2] 

As an example, consider the Werner-like mixture of a W 
state and the maximally mixed state of four modes with a 
single excitation, p mm = ( |0001) (0001 1 + |0010)(0010| + 
|0100)(0100| + |1000)(1000|)/4, 



Pl ( P )=p\W 1 )(W 1 \ + (l-p) Pn 
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Using the above criterion for A(pi(p)) = 3/4 — 3p 2 /4, we 
find that for p > 2/3 we can detect genuine four-mode en- 
tanglement, and for p > v3/3 ~ 0.577 we detect at least 
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FIG. 2: Minimum variances A for various types of four-mode states 
containing exactly one excitation (photon). 



three-mode entanglement. Moreover, for any p < 1, the state 
Pi(p) is entangled, even if just two-mode entangled. 

If the number of modes N is arbitrary, then the minimum of 
A(/9i) for biseparable (N — l)-mode entangled states can be 
shown, after some algebra, to be given by (27V — 3) /N(N — 
1). In the limit of large N, this bound rapidly approaches 
zero, hence making it practically impossible to distinguish in 
this way genuine iV-partite entanglement from mere (N — 1)- 
partite entanglement. 

Finally, we note similar uncertainty-based entanglement 
criteria can in principle be applied to all types of TV -mode 
states with some fixed total number of excitations. If the total 
photon number is larger than 1, however, the unitary transfor- 
mation from product states to a basis consisting of entangled 
states can in general not be performed with linear optical op- 
erations only. Therefore, measurements in such a basis would 
no longer be deterministic in that case. 



III. DETECTING W STATES IN AN EXPERIMENT 

Due to experimental imperfections, an actual state, pro- 
duced in a laboratory, is never a pure state with a fixed num- 
ber of excitations, such as, say, Eq. Q. In experiments with 
atomic ensembles a state pw is routinely generated whose 
single-excitation part has a large overlap with a W state, but 
which contains a significant contribution from the vacuum and 
from states with more than one excitation. As a conserva- 
tive estimate, we can ignore coherent superpositions of states 
with different numbers of excitations [one can get rid of such 
coherences by local operations, see ifFTl l. and hence we can 
write down pw in a generic form 



pw = ppo + qpi + 0--p- q)p>2, 



(ii) 



where the subscript indicates the number of excitations. Typi- 
cally, the magnitude of 1 — p — q is of the order of 1 % or even 
less. The main source of contamination to the desired single- 
excitation part is the vacuum. Moreover, p\ is not necessarily 
a pure state, and is not necessarily a state of N single modes 
either. 

Even if the uncertainty measure from the preceding Section 
would identify the presence of four-mode entanglement in the 
state pi, this does not guarantee that pw itself carries any en- 
tanglement. The standard counterexample |[T8"1 is a four-mode 
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state of the (unnormalized) form 



experiment one would really project onto 4 states of the form 



|+}cx(|0} + e |l)r 4 , 



(12) 



for which the one-excitation part is genuinely four-party en- 
tangled, although the state itself is fully separable. Therefore, 
in order to justify the presence of entanglement in an exper- 
iment it is not sufficient to measure only the variance A(pi) 
of the single-excitation part of the density matrix, but it is cru- 
cial to additionally measure the numbers {p, q}. Once p, q and 
A(p!) are measured one can check if there exists a completely 
separable or biseparable state p tes t with the same values of 
p, g,A(pi). If no such state exists, then one can correctly 
conclude that p\y is entangled. 

More precisely, for fixed values of p and q we wish to find 
the minimum possible value for the variance A, A m j n , con- 
sistent with the various sorts of biseparable or fully separable 
states. In the following we will plot results for the case where 
q = 0.1, which is the relevant case for the experiment lfl3l . 
We will find A m j n in that case as a function of r := 1—p — q. 

Before discussing in turn the various classes of separable 
and biseparable states, we make several remarks: 

We note that A m i n within each such class cannot increase 
with decreasing q. The reason is that given any state p we can 
always mix in the vacuum po without changing the variance 
A, and without increasing the entanglement. But this mix- 
ing operation clearly does decrease q. Hence A m i n cannot 
increase with decreasing q. 

Similarly, we could mix in a fully separable state containing 
more than a single photon in some given mode, e.g., a tensor 
product of the vacuum and one mode with 2 or more photons. 
This again does not affect A, and does not increase entangle- 
ment, but does decrease q. For this reason, in our attempts to 
find the minimum variance we do not have to consider states 
with more than a single photon in any given mode, as those 
states will have a larger value of A than the minimum possi- 
ble for given q. 

Moreover, we could take a state of N single modes and 
convert it into a state of multiple modes in N locations, by 
locally applying a random unitary operation. This local oper- 
ation does not move a state up the entanglement hierarchy and 
does not affect any of the quantities A, q and r. Thus, exclud- 
ing fully separable states and biseparable states of N single 
modes is sufficient for detecting entanglement. 

Because p 1 is subnormalized to q, we have, instead of the 
inequality <|3j, the inequality 



gA(pi) > p m g m A(jO m ,i), 
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where q = ^2 m p m q m and p = J2 m PmP m - 

Finally, instead of projecting onto the 4 states (|4|l-((7]), in an 
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and vary over all three phases <f>k,k = 1,2,3 (which one 
accomplishes by inserting phase shifters in the appropriate 
modes) to find the minimum variance, thus optimizing the 
entanglement test. Our method is otherwise independent of 
which values of <f>k attain that minimum. 



A. Fully separable four-mode states 

It is relatively easy to account for all separable and bisep- 
arable states in the case of four modes. Let us first calculate 
p, q, A(pi) for fully separable states. We first consider pure 
states 1 ip s ) of the form 



A jjejl)) 



(15) 



for complex parameters e^. As argued above, we do not have 
to consider states with more than a single excitation in any one 
mode. For the pure state |^ s ) a corresponding density matrix 
can be constructed from \ip s )(i/j s | : 



Ps = PPo + QPi + rp>2, 



(16) 



where 



and 



4 = 1 

We can visualize the set of pure completely separable states, 
and in particular its border, by plotting values of A(p s l ) ver- 
sus r for a fixed value of the single-excitation probability q, 
by randomly varying over all values of e, consistent with that 
value of q. By symmetry, it is clear the minimum variance 
will be obtained for real parameters. The result is shown in 
Figure [3] and we can clearly identify the region of full separa- 
bility, the lightly shaded area (colored in yellow, online). The 
minimum value of A(p s i) at r = is 3/4 which is an agree- 
ment with our previous discussion (see Fig.|2]i. It is instructive 
to point out once again that, even though A(/? S) i) approaches 
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FIG. 3: (Color online) Scatter plot (in yellow) of the variance 
A(p Si i) of the single-photon part for randomly chosen, pure, fully 
separable states versus the probability of finding multiple excitations, 
r for a fixed single-photon probability q = 0.1. The black crosses 
are data points for a particular subset of pure states, attaining the ex- 
tremum values of the variance for the set of pure states. Also plotted 
is the variance for randomly chosen mixed states (in green). For this 
particular value of q those values for the variance fall within the con- 
vex hull of the graph for pure states (the red line is the convex hull 
of the black curve). The region below the lowest red line then corre- 
sponds to entangled states (as indicated by the word "ENTANGLE- 
MENT"): but this includes two-mode, three-mode and four-mode 
entanglement. 



zero for sufficiently large values of r [see also the last subsec- 
tion of this Section], the density matrix p s is and remains fully 
separable [cf. the example mentioned above]. 

Moreover, we managed to find the pure states living on 
the pure-state boundary, indicated in black in Figure [3] The 
boundary can be parameterized by two parameters, either q 
and r, or, more simply, by e and e. Namely, the extremum 
values of the variance for pure fully separable states turns out 
to be attained for states of the form 

|</V>«(|0) + £ |l))(|0)+e|l)r 3 . 

Now one notices the lower border for pure states in Figure [3] 
is not convex as plotted. This may indicate that points corre- 
sponding to certain mixed states may fall below the pure-state 
boundary. Thus we have also tested randomly chosen mix- 
tures of random pure states as well as mixtures of states on 
that boundary. And some mixed states (plotted in green) in- 
deed have a smaller variance. Thus, the minimum variance 
is attained by mixed states in this case, and the correct lower 
bound is indicated in red. This lower bound coincides with 
the convex hull of the graph for pure states. 



B. Biseparable states with at most two-mode entanglement 

The next class of states to consider is biseparable states 
i.e. states that can be described by a density matrix p\>i S — 



M 

Tli pt®pf ■ The division into subsystems A and B in the case 

i=l 

of four modes has two distinct possibilities - either system A 
represents one of the modes and system B consists of the re- 
maining three modes or both systems A and B represent two 
modes each. We will study the latter case first. We represent 
a pure biseparable state with at most two-mode entanglement 
by 

\iI>)ab = W)a®\$)b, (17) 

with both two-mode states \tp)k for k — A, B of the form 

|^> fc cx|00)+e fc |01) + 4|10), (18) 

where we now included phase factors into the parameters t k 
and e' k . For the same reason as given in the preceding subsec- 
tion, we do not have to consider contributions from terms with 
more than a single excitation in any one mode. The expres- 
sion for the variance A is symmetric under the interchange of 
any two modes, and so it is immaterial which two modes con- 
stitute system A. We again vary over the complex parameters 
efc, e' k for fixed value of the single-photon probability q = 0.1 
to find the set of all pure biseparable states with at most two- 
mode entanglement, as a function of the multiple-excitation 
probability r. By symmetry the minimum variance is attained 
for real coefficients. 

The lightly-shaded (yellow) area in Figure |4] then depicts 
the set containing all biseparable states with at most two-mode 
entanglement. Indeed, we have checked explicitly that points 
corresponding to mixed states fall within the shaded region, 
unlike in the preceding case of fully separable states. The 
shaded region of Figure |4] includes that of Figure [3] simply 
because the set of fully separable states is a subset of the set 
of states with at most 2-mode entanglement. The minimum 
value of A at r — is 1/2, confirming the result from Fig. [2] 
Just as in the preceding subsection we find the pure states liv- 
ing on the boundary. The boundary is again parameterized by 
two parameters, e and e. Namely, the minimum variance is 
attained for biseparable states of the form 

\A,i)AB cc (|00) + e|10) + e|01))(|00) + e|10) + e|01». 

In this case it is straightforward to extract the minimum vari- 
ance as a function of q and r: 

A min = 1/2 - 2r(l - q)/q 2 + 2r 2 /q 2 , (19) 

which is indeed almost linear in r when r<Cl. Moreover, this 
boundary is convex. We can rewrite the minimum variance 
more compactly as 

A min = l/2-2rp/q 2 . (20) 

C. Biseparable states with three-mode entanglement 

A pure state of the entire four-mode system (with up to two 
excitations) that has at most three-mode entanglement can be 
described by the following biseparable vector, 

|^>AB CX (|0>+€i|l» ® 

(|000) + e a |100) + e 3 |010) + e 4 |001)), (21) 
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FIG. 4: (Color online.) Same as Figure [3] but for pure biseparable 
states with at most two-mode entanglement. The graph is convex, 
and points corresponding to mixed states (plotted in green) fall within 
the yellow region. The region below the black curve corresponds to 
at least three-mode entanglement. 



FIG. 5: (Color online.) Same as Figure|3](g = 0.1), but for bisepara- 
ble states with at most three-mode entanglement. Points correspond- 
ing to mixed states fall within the lightly shaded (yellow) region, and 
are plotted in green. The region below the black curve corresponds 
to states with genuine four-mode entanglement. 



where we have arbitrarily chosen the first mode to be the sys- 
tem A. The analysis, however, is symmetric with respect to 
our choice for the system A. In the second term we do not 
have to consider states with more than a single photon in sys- 
tem B. Although the measurement determining whether there 
are multiple excitations in the three modes comprising system 
B is not a local filtering operation in the usual sense, it is local 
with respect to the bipartite cut A vs B, which is the relevant 
cut in this case. 

The result for q = 0.1 is plotted in Fig. [5] Like in the case of 
separable states we observe the existence of points where the 
values of A (pi) are close to zero. Again we should not mis- 
judge the presence of entanglement in these states, since the 
states we are operating with are biseparable by construction. 
The region for r larger than w 2 x 10~ 3 where the minimum 
variance no longer is a decreasing function of 7-, contains no 
physical states with smaller variance. 

The lightly-shaded (yellow) region depicts the convex set of 
biseparable states with at most three-mode entanglement, and 
includes the set of fully separable states, although not nec- 
essarily the set of states with two two-mode entangled states 
[it does for q = 0.1]. We have explicitly verified that points 
corresponding to mixed states (plotted in green) fall within 
the yellow region. The minimum value of A at r — is 
perhaps a little hard to discern, but is indeed equal to 5/12, 
the value obtained analytically in the preceding Section. The 
lower boundary (plotted in black) corresponds to states of the 
form 

|V>>Ai?oc(|0) + e|l))<g) 
(|000) + e|100) + ejOlO) + e|001», (22) 

with real and positive e > e. Since the boundary of minimum 
variance is the lowest for this type of biseparable states, it is 
the relevant boundary for the purpose of detecting genuine 4- 
mode entanglement. For this reason we plot these boundaries 
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FIG. 6: (Color online.) Same as Fig.|5]but for q = 0.4 

for several values of q. For increasing values of q the mini- 
mum possible variance for 3-mode entangled states increases 
and reaches the limit of max(A m j n ) = 5/12 for q — ► 1. 
Figures |^|8] approach this limit for values q = 0.4 through 
q = 0.7 to q = 0.9. 



D. Putting it all together 

Based on an exclusion analysis, a practical inseparability 
criterion can be formulated. In an experiment aimed at de- 
tecting a genuinely four-mode entanglement one measures the 
values of r, q and A. Then one plots, according to the previ- 
ous considerations, values of A versus r for all separable and 
biseparable-state models, feeding in the value of q attained 
from the experiment. The measured values of p, q and A are 
represented by a single point in that plot. If that point lie out- 
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FIG. 7: (Color online.) Same as Fig.|5]butforg = 0.7 
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FIG. 8: (Color online.) Same as Fig.|3]but for q = 0.9. 



side all three shaded regions of the model plots, the state pro- 
duced in the experiment must carry genuine four-mode en- 
tanglement. Partial conclusions about entanglement can be 
reached when the point falls outside some and inside other 
regions. In particular, if the measurement point lies outside 
the shaded region of Fig. [3] but inside the shaded regions of 
Figures [4] and [5] one can only conclude one has an entangled 
state, but it could be merely two-mode entangled. If the point 
falls outside the shaded regions of both Figures [3] and |4] but 
inside the shaded region of Figure [5] one has at least three- 
mode entanglement. Of course, if the point falls inside the 
shaded region of Figure [3] no firm conclusion can be reached 
about entanglement, as there is a fully separable state consis- 
tent with one's values for p, q, A. 

We plot the three minimum-variance boundaries for dif- 
ferent small values of q as a function of a scaled variable 
R := 8rp/3q 2 (see the next subsection for an explanation 
for this choice of variable) in Fig. [9] One sees the boundaries 
depend only weakly on that parameter. 




0.2 0.4 0.6 0.8 1 

R 



FIG. 9: (Color online.) Boundaries for the minimum variance for the 
three types of biseparable states as functions of R := 8rp/3q 2 for 
10 values of q = 0.02, 0.04, . . . 0.2. The reason for choosing this 
particular variable is given in Section [TlI E| The lowest-lying (green) 
curves correspond to 3-party entangled states, the highest-lying (red) 
curves correspond to fully separable states, the middle (blue) curves 
correspond to biseparable states with 2 -mode entanglement. The 
variance depends only weakly on q for the red curves, and is inde- 
pendent of q for the green and blue curves. 

E. Some necessary conditions for entanglement 

Let us finally consider the conditions on entanglement in 
the simple situation where the variance A(pi) vanishes and 
where q is not too large. We consider the same three classes 
of unentangled states as before. 

1. Fully separable states with A(pi) = must be of the 
form (|0) + e|l))® 4 . For such states, the point at which 
the variance is zero is characterized by 

8 rp q q 2 

3 q 2 6p 96p 

For small values of q, we can give the approximate re- 
lation, valid for fully separable states: 

R>1, 

with R = 8rp/3q 2 being the quantity appearing on the 
lhs of ( |23"| l. A necessary (although not sufficient) con- 
dition for any type of entanglement is then simply 

R < 1. 

For Figure [3] in which we took q = 0.1, this places a 
strict upper limit on r of r < 4.125 x 10 -3 for entan- 
glement to be detectable through A. 

2. Biseparable states with A(pi) = and at most two- 
mode entanglement must be of the form (|00) +e|01) + 
e|10))® 2 . For such states, the boundary of zero variance 
is at R = 2/3, and hence all biseparable states satisfy 

«4 
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For Figure |4j in which q = 0.1, this places a strict up- 
per limit on r of r < 2.75 x 1CP 3 for entanglement 
involving at least three modes to be detectable through 
A. 

3. Biseparable states with A(pi) = and at most three- 
mode entanglement must be of the form (|0) + e 1 1} ) (g> 
(|000) +e|001) + e|010) +e|100). For such states, we 
similarly derive 

For Figure [5] in which q = 0.1, this places an upper 
limit on r of r < 2.06 x 10~ 3 for entanglement to be 
detectable through A. 

In order to demonstrate genuine four-mode entanglement one 
must violate all of these conditions. That is, one must violate 
the strongest of these conditions, and hence one must have 



R < 



1 



(24) 



This condition for four-mode entanglement is necessary but 
not sufficient for nonzero values of A(pi). The form of the 
conditions also indicates why the scaled variable R, used in 
Fig. [9] is a useful quantity for small q for fully separable states, 
and for biseparable states irrespective of the value of q. 



IV. LOSSES AND ASYMMETRIES 

So far we have assumed that the variance measurement de- 
vice is ideal: beamsplitters (see Fig. [T} were assumed lossless 
and perfectly balanced, and detectors were perfect. In this 
subsection we relax those conditions and describe the modi- 
fications necessary to include these imperfections. First, we 
consider the effect of imbalanced beamsplitters. 



A. Imbalanced beamsplitters 

Suppose, then, we have the same setup as depicted in Fig. 
[T| but with the 4 beamsplitters having reflection and transmis- 
sion probabilities \tk\ 2 and \rk\ 2 not necessarily equal to 1/2. 
Consider one output mode, say the top one. There is one path 
a photon can take from input mode 1 to reach the top output 
mode: it has to reflect off of beamsplitter 1 and it has to re- 
flect off of beamsplitter 3. The amplitude for that path is then 
r\rz in terms of the reflection amplitudes of beamsplitters 1 
and 3. Here we ignore phase factors due to propagation (they 
can be trivially inserted in the end). Similarly, a photon from 
input mode 2 can reach the top output mode along just one 
path, with amplitude iir 3 . Writing down the amplitudes for 
photons starting in input modes 3 and 4 shows that a photo- 
detection at the top output mode projects onto the (input) state 



This is a properly normalized state, even if the beamsplitters 
are not balanced. The normalization follows from the relation 
\ r k\ 2 + \tk\ 2 = 1 for lossless beamsplitters. 

We can similarly write down the states onto which one 
projects if detecting a photon in one of the remaining output 
modes: 

\W 2 ) = £ir 4 |1000) +rir 4 |0100) + i 2 i 4 |0010) 

+r 2 i 4 |0001), (26) 

\W 3 ) = r^llOOO) +i 1 t 3 |0100) +r 2 r 3 |0010) + 

£ 2 r 3 |0001>, (27) 

\Wi) = £ii 4 |1000> + r!t 4 |0100) + i 2 r 4 |0010> 

+r 2 r 4 |0001). (28) 

These states, too, are normalized. Moreover, the 4 states 
are all orthogonal, as follows from the (unitarity) relation 
t* k rk + tkr\ = 0. One can still calculate the variance of pho- 
todetector counts, using the modified projectors onto the W 
states, but that variance will not give as much information as in 
the balanced case about four-mode entanglement. For exam- 
ple, consider the extreme case of a mirror replacing beamsplit- 
ter 4: that is, assume now that r 4 = 1 and i 4 = 0. Then states 
IW2) and IW4) are no longer four- mode entangled states, but 
only two-mode entangled states. Thus, certain two-mode en- 
tangle states would give rise to a zero variance in this extreme 
case. 

This implies that even if one's experiment cannot use per- 
fect 50/50 beamsplitters, one should at least try to make them 
as balanced as possible. In such cases one needs in general a 
lower variance A than in the ideal balanced case to conclude 
one has four-mode entanglement. 



B. Losses 

Now consider losses. We can model linear losses (both 
propagation losses, and inefficiencies of the photodetectors) 
by imagining lossless paths but with additional beamsplitters 
reflecting away some portion of the light in the lossy paths. 
The output of those additional beamsplitters does not lead to 
the output detectors, but to other (unmonitored) output modes. 
The overall transformation from input to output is still unitary, 
which implies there must also be additional input modes (just 
as many as there are unmonitored output modes). A photode- 
tection in one of the desired output modes projects onto a set 
of orthonormal states on the larger Hilbert space of all input 
modes. If we write down the projections of those states onto 
the 4 input modes of interest, we will end up with subnormal- 
ized states. For example, considering for the moment (see the 
next subsection where we take into account multiple excita- 
tions) only states with exactly one photon, a detection in the 
top output mode projects onto the state 



\W{) = TunrallOOO) 
+T 31 r 2 i 3 |0010) 



■ r 2 itir 8 |0100) 
-T 41 i 2 i 3 |0001) 



(29) 



\Wi) =r 1 r 3 |1000) + t 1 r 3 |0100) 



r 2 i 3 |0010) + i 2 i 3 |0001). 

(25) 



where the transmission amplitude Tf-i is the product of all loss 
amplitudes encountered by a photon propagating from input k 
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to the top output detector (this includes the inefficiency of the 
detector). 

The variance we are interested in is conditioned on detect- 
ing (at least) one photon in the desired output modes. Once 
we detect a photon in the top mode, we do renormalize the 
state \W[) and project onto: 



\W{') = 



\W{) 



(W[\W{) 



(30) 



The 4 states onto which we project conditionally, \W[!) for 
k = 1 ... 4 are, therefore, properly normalized, but they are 
not orthogonal, unless all losses are balanced, i.e., if = 
constant for all/, k = 1 ... 4. 

Again, we still can use a variance based on the modified 
nonorthogonal projectors, but that variance will give less in- 
formation than in the ideal lossless balanced case. For in- 
stance, if all photodetectors but one are completely inefficient 
and never detect any photon, the variance would be zero for 
any input state. Thus, in an actual experiment one would have 
to make the losses as balanced as possible in order for the 
variance to contain as much information about four-mode en- 
tanglement as possible. Of course, one would also like to limit 
the size of the losses for various different reasons. 

With the new projectors onto the nonorthogonal states 
\Wf!) in hand, we can perform the same calculations as we 
did in the ideal case: find the minimum variance consistent 
with unentangled input states, input states with two-mode en- 
tanglement, and input states with three-mode entanglement, 
respectively, for fixed values of q and r. We display three 
illustrative examples (for q — 0.1): in Fig. [T0| we assume 
no losses but unbalanced beamsplitters, in Fig. 1 1 we assume 
losses, but balanced beamsplitters, and in Fig. 12 we show 
the effects of the combination of both losses and imbalances. 
All figures show the tendency of the minimum variance to de- 
crease compared to the ideal lossless and balanced case. 



C. Measured variance vs A (pi) 

In the presence of losses the measured variance, A m , is 
not just due to the single-excitation part, but from the multi- 
excitation part of the input state as well. Thus, the measured 
variance has to be corrected (upwards in fact) in order to find 
an estimate for the variance A(pi) due to the single-excitation 
part, because that is the quantity we used above to detect en- 
tanglement. 

We discuss a simple case (balanced losses throughout the 
system and the use of non-number resolving threshold detec- 
tors l22ll ) where we find we simply have to multiply the mea- 
sured variance with a factor c > 1 to obtain an upper bound 
on A. That is, the variance A is upper-bounded by cA m . We 
now evaluate c. 

Consider the propagation of the purported experimental 
state (here we change notation to make it easier to keep track 
of the meaning of all symbols) 
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FIG. 10: (Color online.) Minimum variance curves for the case 
where all beamsplitters are 55/45 (dashed lines) rather than 50/50 
(solid lines). 
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FIG. 1 1 : (Color online.) Minimum variance curves for the case 
where there is one lossy path with transmission probability of 60% 
(a typical parameter) (dashed lines), compared to the ideal lossless 
case (solid lines). 



Under balanced losses (which can be characterized by a single 
transmission efficiency \T\ 2 ) this state transforms to pt where 




PT =PoPo +p' 1 {q' 1 p i i ) + (1 



9i)Pi 2) ) 



■p' >2 p { *\. (32) 



ill (r) 

Here, p\ is the 1 -photon subspace of p\ T <2 originating from 



the z-photon pi subspace of pZl for i £ {1,2}, and /?> I is 

the 2-photon subspace after the transmission. To the leading 

(2) 

order of p^ (neglecting 3-photon and 4-photon subspaces), 
are 



Pw = PoPq + P1P1 + P>2P>2- 



(31) 



Po=Po + (lHH 2 ) P i + (l-|T 



2\2 



P2 



(33) 
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FIG. 12: (Color online.) Minimum variance curves for the case 
where there is one lossy path with a transmission probability of 60% 
and unbalanced 55/45 beamsplitters (dashed lines), compared to the 
ideal lossless and balanced case 



pi = |T| 2 Pl + 2|T| 2 (1 - |7f)p 2 



i>2 



and q[ is given as 



<h 



in 



in 



P2, 



Pi 



|T|2 pi + 2|T| 2 (l-|T|2)p 2 - 



(34) 



(35) 



(36) 



Thus, if we denote by P& the normalized probability of 
the detector (assumed to be non-number resolving) in mode 
k finding (at least) one photon , then we have 



Pk = 



Mip(i) 
Q l ' k 



Pi(l 



9i) p (i) 

Q 2 - k 



P2 p (2) 

Q 2 ' k ' 



(37) 



where Q — p[ + p' 2 . Here, PP") (Pi V) is the probability of 
a 1 -photon in output mode k originating from the l(2)-photon 

(2) 

subspace p\ (pi), and P 2 k i s me probability of 2-photon in 
output mode k. To be conservative (for our purposes of finding 
a sufficient condition for entanglement), we assume that the 
two photons are directed towards one detector at a time so 
that we cannot distinguish P} 1 / from P?l. By denoting qy = 
pWi/Q as tne probability of detecting desired events and Xk 
as the normalized probability of detecting undesired events 
(that is, P 2 (1 2 and P^X we get 



Pk = qiP§ + (1 - qi)X k . 
The measured variance A m is given as 

A m = l -J2 Pi 



(38) 



(39) 



On the other hand, the 1 -photon variance A is defined as 

(1)\2 



a(pi) = i-EO 



(40) 



As a conservative correction to A m , we assume that the un- 
wanted events (Xf.) are all directed towards the output mode j 
which contains the maximum 1-photon probability P^J (i.e. 
Xj = 1 and Xk — for k ^ j). This way, the measured 
variance is lower than the variance A. Thus our conservative 
bound gives then 

A m = 1 - ( 9l P« + (1 - qi )f q \ Y^(P^) 2 - (41) 



1 Using the inequality 2(1 - P^) > (1 + P x ( ^ )(1 - P$) 



1 - (P^) 2 , we obtain 



A m > 9 i[A + (1 - 9l ) ^(P^) 2 ] > qiA. 
Therefore, we obtain a correction factor of 



1 



where 



9i = 



Pi 



(42) 



(43) 



(44) 



(45) 



Q p! + (2-|T| 2 )p> 2 ' 
In the limit of po w 1, the correction factor becomes 

c«l + ^(2-|T|VP. 
V. SUMMARY AND DISCUSSION 



We demonstrated how to verify A^-party entanglement of W 
states or states lying close to W states, in the case quantum in- 
formation is encoded in the number of excitations per mode. 
Our method takes into account the presence of the vacuum 
state, as well as multiple excitations; moreover, it takes into 
account losses during the verification measurements, as well 
as imperfect beamsplitters. The method applies to any number 
of modes, but we focused on four modes for illustrative pur- 
poses, as the method was applied in an actual experiment lfT3l 
to four modes. A relatively straightforward set of measure- 
ments allows one, in that case, to distinguish genuine four- 
party entanglement from three-party entanglement, which in 
turn can be distinguished from two-party entanglement and 
fully separable states. One must obtain estimates of three pa- 
rameters: a variance A determined by the single-photon part 
of the state, the single-photon probability q, as well as the 
multi-photon probability r. For example, the simple condition 



of Eq. ( 24 1 is a necessary condition for genuine four-party en- 
tanglement (where our definition of genuine multi-partite en- 
tanglement is more severe than usual) which involves only r 
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and q. To obtain sufficient conditions one must also include 
the value of A in the analysis. 

The measurement of A combines the various modes by 
simple beamsplitters, and is thus a nonlocal measurement. In 
this way one does not need local oscillators, which one would 
need if the entanglement verification method used local mea- 
surements only [19|. In our case the modes interfere with 
each other, rather than with external reference beams. Thus, 
our method cannot be applied to eliminate local hidden vari- 
able models (through Bell inequalities, for example, in the bi- 



partite case), but it can be applied to verifying entanglement, 
which is a very different beast indeed 0. 
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